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SEEDING THE FORMATION OF COLD GASEOUS CLOUDS IN MILKY WAY SIZE HALOS 



OS 

o 
o 

(N 



O 

u 

6 

(N 
> 

00 

wn 
O 

o> 
o 



X 



Dusan Keres & Lars Hernquist 

Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138; dkeres9cfa.harvard.edu 

Draft version July 1, 2009 

ABSTRACT 

We use one of the highest resolution cosmological SPH simulations to date to demonstrate that cold 
gaseous clouds form around Milky Way size galaxies. We further explore mechanisms responsible for 
their formation and show that a large fraction of clouds originate as a consequence of late-time fila- 
mentary "cold mode" accretion. Here, filaments that are still colder and denser than the surrounding 
halo gas, are not able to connect directly to galaxies, as they do at high redshift, but are instead 
susceptible to the combined action of cooling and Rayleigh- Taylor instabilities at intermediate radii 
within the halo leading to the production of cold, dense pressure-confined clouds, without an associ- 
ated dark matter component. This process is aided through the compression of the incoming filaments 
by the hot halo gas and expanding shocks during the halo buildup. Our mechanism directly seeds 
clouds from gas with substantial local overdensity, unlike in previous models, and provides a channel 
for the origin of cloud complexes. These clouds can later "rain" onto galaxies, delivering fuel for star 
formation. Owing to the relatively large cross section of filaments and the net angular momentum 
carried by the gas, the clouds will be distributed in a modestly flattened region around a galaxy. 
Subject headings: galaxies: formation — Galaxy: formation — Galaxy: halo — instabilities — meth- 
ods: numerical 



1. INTRODUCTION 

Observations of neutral hydrogen around the Milky 
Way (MW) have revealed that our Galaxy is surrounded 
by a large number of cold gaseous clouds. Owing to 
their high velocities compared to galactic rotation, it 
was proposed that these objects reside outside the disk 
(jMuller et al.l [19631 . and they became known as high- 
velocity clouds (HVCs). The masses of most HVCs 
are uncertain because of their uncertain distances (sev- 
eral HVC complexe s are now placed a t D ~ 10 kpc 
(|Wakker et al.l l2008t iThom et all l200l)), but are es- 
timated to be 10 5 - 1O 7 M (|Wakker fc van Woerdenl 
1997). Similar HI clouds and cloud complexes have 
been found aroun d other galaxies (for a review see 
iSancisi et aLll2008[ ). 

The overall kinematics of HVCs implies net inflow, 
motivating a picture wher e HVCs rep resent the build- 
ing blocks of our Galaxy (jOortJ [l970n and provide the 
supply of gas needed for its relatively high star forma- 
tion rate. In fact, our Galaxy has been forming stars 
at a steady rate for the last several Gyrs (|Binnev et all 
2000), implying that the gas consumed must be contin- 
uously replenished. 

Many HVCs can be ionized by the external UV field, 
so HI observations likely detect only the densest of these 
structures that are able to self-shield. The presence of 
clouds around other galaxies is also implied by absorp- 
tion line studies at both l ow and intermediate r edshifts 
(jLanzetta fc Bowed 119901: IGauthier et~aTI 12009ft . which 
typically find gaseous clouds to distances of several tens 
of kiloparsecs in a halo. It is likely that these ab- 
sorbing clouds are analo gs to the local HVC population 
(jMaller fc Bullock! l200l . 

Despite numerous theories, the origin of these clouds 
remains uncertain. Proposed mechanisms fall into the 
following three classes: 1: Cooling instabilities in the hot 



medium (jFieldl Il965t iMo fc Miralda-Escudd Tl 996ft that 
can interact with the hot halo gas resulting in net infall 
onto the Milky Way (jMaller fc Bullockl 12004) . 2: Pro- 
duction by a "galactic fountain" where outflows driven 
by supernovae from the disk c ondense in the h ot corona 
and fall back onto the MW (iBregmanl 11980ft. possibly 
entra ining hot gas from the halo (|Fraternali fc Binnevl 
2008). 3: Interaction of hot h alo gas with infalling satel- 
lites (e.g. lPutman et aLll2003ft . to account for clouds with 
specific properties; e.g. those in the Magellanic Cloud 
system. In principle, all of these channels could operate, 
complicating modeling and data interpretation. Here, we 
concentrate on cloud formation from cooling and hydro- 
dynamical instabilities in infalling gas. 

Cosmological hydrodynamic simulations have shown 
that at high redshift the bulk of the gas supplied to galax- 
ies is accreted from the int ergalactic medium through fil- 
aments in a "cold mode" ([Katz et all 120031: iKeres et al l 
20051: iKeresI 120071: jKeres et all 120091: Dekel fc Birnboiml 
200llOcvirk et al ll2008l : iBrooks et al.ll2009ft even after a 



hot halo atmosphere develops. Owing to lower densities, 
at late times there is a shift in the physics that gov- 
erns the accretion, where filamentary flows get disrupted 
in the central parts of massive halos. We examine this 
late-time accretion and show that these disrupted fila- 
ments are still able to supply cold gas to MW size galax- 
ies. We demonstrate that accretion through filaments 
creates density inversions in halos, triggering a combina- 
tion of cooling and Rayleigh- Taylor (RT hereafter) insta- 
bilities that produce clouds, supplying MW like galax- 
ies with gas. Earlier work hinted that the "leftovers" 
from filamentary might contribute to clou d formation 
(iSommer-Larseni 120061: jConnors et al.l 120061: iPeek et all 
l2008 UKeres et al.l 12009ft but did not demonstrate why 
and how this process operates. 



2. SIMULATIONS 
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We use the Gadge t- 2 Smoothed Particle Hydrody- 
namics (SPH) code (jSpringell [2005) which simultane- 
ously conserves energy and entropy, preventing numeri- 
cal phase mixing. This is potentially an important prop- 
erty for modeling cloud formation as many SPH implc- 
mentations can mix cold, dense and ho t, dilute phases 
(jSpringel fe Hernquistl 120021 : lKeres!l2007h . Gravitational 
forces are calculated combining Particle Mesh and hi- 
erarchical tree algorithms. In order to capture the for- 
mation of ~ 10 6 M Q clouds at high spatial resolution in 
a large volume, we adopt a zoom-in technique, where 
we first simulate a cosmological box of lO/i -1 Mpc on 
a side using only dark matter (DM hereafter) particles. 
Next, we select a region that will collapse into a MW 
size object by z=0. We populate this region and its 
surroundings at high resolution with gas and DM par- 
ticles, introducing ACDM density fluctuations down to 
the new inter-particle separation scale. Outside of this 
region, the density field is sampled more coarsely with 
high mass DM particles to provide long range forces. 
Suc h initial conditio n s are re-generated at z = 99 using 
the lEisenstein fc Hul (|1999f ) power spectrum in a ACDM 
cosmology: fi b = 0.044, Q, m = 0.27, h = 0.7, a s = 0.8. 
We include standard cooling processe s for a gas of pri- 
mordial composition (iKatz et al.lll996h and the UV back- 
ground (jHaardt fc Madaulll996[ K For simplicity, we ne- 
glect galactic outflows. 

In total, the simulation has 6.6 million particles. The 
most massive halo in the high-resolution region (here- 
after, MWH), analyzed in detail in what follows, at z=0 
has Mhaio = 7 x lO n M0 and contains 1.4 million bary- 
onic and 0.9 million DM particles. Each gas particle 
begins with a mass of 9 x 10 4 h~ 1 MQ, while stars are half 
this mass. Gravitational softening is the same for DM 
and baryonic particles, e — 400ft~ 1 (l + z) _1 pc (Plum- 
mer equivalent). 

3. RESULTS 
3.1. Clouds 

In Figure [TJ we show the inner regions of MWH. Gas 
particles with nn > 4 x 10~ 4 cto~ 3 are plotted close to 
face-on and edge-on projections of the disk, revealing the 
presence of a population of dense gaseous clouds together 
with a central galactic disk. Plotting gas colder than 
30000K selects an identical cloud population. At low red- 
shift MWH contains a large number of cold, T c ~ lO 4 !^, 
dense clouds in a 40/i _1 kpc (edge on) times 60ft, -1 kpc 
(face on) region around the disk. We plot gas particles 
with small arrows to indicate the velocity field, demon- 
strating that the clouds exhibit a global sense of rotation, 
similar in orientation to the outer disk. 

We identify clouds with T < 30000 K and n H > 0.004 
cm~ 3 using the code SKID. The results are mainly de- 
termined by the temperature threshold and are robust 
for any T < 10 5 K. To distinguish clouds from galaxies 
we modified SKID to include thermal support of gas to 
properly identify self-gravitating systems. Bound objects 
include the main galactic disk, several satellite galaxies 
(one visible in the upper left corner of Figure [T]) , and a 
small number of gaseous clumps that are not associated 
with a DM or stellar component. These self-gravitating 
clouds are close to the halo center, where the external 
pressure compresses them to high densities. For sim- 



plicity, we exclude all bound objects from the remainder 
of the analysis, but apart from the main gaseous disk 
they contribute negligibly to the cold gas budget. Cloud 
masses range from 10 6 to few 10 7 M©. We note, how- 
ever, that we arc limited at low masses by the simula- 
tion resolution. Some of the most massive clouds are 
at larger radii and often represent a dense slab of gas, 
which might be susceptible to further instabilities. Typ- 
ical column densities (for the total cold gas, i.e. no ion- 
ization correction which could be substantial) of clouds 
are JV ff ~5x 10 18 - l^cmT 2 . 

3.2. Gas in the inner halo 

The amount of virialized gas increases with time as 
the halo grows and hot gas fills l ow density region s in- 
between colder, denser filaments ([Keres et~all[2005l ). At 
high redshift, filaments survive within hot halos even 
when the hot medium dominates the halo gas budget 
(M h > 3 x 10 n M Q ). However, at lower redshift the den- 
sity of filaments decreases and the temperature increases 
so they are more easily destroyed at small radii by large 
scale shocks in a hot halo. This happens where halo 
shock conditions are favorable for hot, overd ense gas to 
propagate outwards (jBirnboim fc Dckel 2003) . Filamen- 
tary remnants that penetrate into halos but not quite to 
galaxies set the stage for cloud formation. For the MWH 
this epoch starts roughly at z = 1 — 1.5 after the last 
major merger is completed. The last significant (~ 1 : 7) 
merger occurs around z=0.5 and soon afterwards MWH 
reaches a quasi-static configuration that lasts 2-3 Gyrs, 
although cloud formation proceeds in a similar manner 
for z < 1. At z ~ around 10 percent of the halo gas (ex- 
cluding gas in the galaxies) within the inner 100/i -1 kpc, 
1.5 — 2 x 10 9 M Q , is in cold gaseous clouds as a result 
of a balance between cloud creation, accretion and de- 
struction. This is consistent with estimates of the total 
mass of the HVC population (|Putmanl 12006). The den- 
sity of the hot component at a distance of 50 kpc from 
the main galaxy is nn ~ .5 — 1.5 x 10~ 4 cm~ 3 , consis tent 
with observational limits (jGrcevich fc Putmanll2009[ ). 

3.3. The origin of the cold clouds 

To understand their origin we follow cloud particles 
backwards in time and determine the maximum temper- 
ature that they had in the past, T max . This tells us if 
the clouds form through instabilities in a hot corona or 
if they were seeded by colder structures. Figure [5] shows 
the results of this analysis (for "resolved" clouds, > 32 
particles) . Most of the mass in the clouds never reached 
0.5 - 1.5 x 10 6 K, characteristic of the MWH at z=0, but 
instead clouds form from, 1 — 1.5 x 10 5 -fT material, before 
they cool to ~ 10 4 K. These are temperatures character- 
istic of low-redshift, large-scale filaments fueling gas onto 
MWH. The distribution of T max is similar in the redshift 
interval < z < 1, implying a connection between cloud 
formation and "cold mode" accretion active at higher 
redshifts. 

We illustrate this process in Figure [3) An overdense, 
warm filament penetrates to intermediate radii in the 
halo, visible in the panels which show density and tem- 
perature. The filament is not directly connected to the 
galaxy but manages to reach 30-40 kpc radius before be- 
ing heated to higher temperatures. The large impact 
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Fig. 1. — Distribution of dense clouds in a 200fo 1 kpc box centered on the most massive galaxy in MWH. Particles are color coded by 
density from riu = 4 X 10~ 4 cm" 3 (darkest blue) to 4 cm -3 (bright yellow). Clouds have temperatures T c ~ 10 4 K. 

so it is important to estimate the timescales for a dom- 
inant mode to affect and destroy such penetrating fil- 
aments before they become cooling or RT unstable at 
small radii. The relevant timescale is: 
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Fig. 2. — Histogram of the maximum temperature cloud particles 
reached before condensing into clouds, in bins of 0.2 log T max . 

parameter with respect to the central galaxy creates a 
density inversion where "heavier" fluid lies on top of a 
"lighter" one in an external gravitational field (of the 
parent halo), providing the conditions for RT instabil- 
ities to operate. The initial instabilities that act on 
overdense filamentary gas often produce cometary tail- 
like morphologies, similar to RT instabilities in the ISM 
(e.g. IP 'dell fc Burkertl fl997h . These are further com- 
pressed by higher pressure gas at slightly smaller radii 
and quickly become cooling unstable, forming cold, dense 
clouds. Additional sites of cloud formation throughout 
the MWH can be seen in the figure, often where com- 
pression of infalling gas by halo shocks and pressure from 
surrounding hot gas promotes faster development of in- 
stabilities. 

The filament is also moving with respect to the sur- 
rounding lower density gas. It is therefore possible that 
Kelvin Helmholtz (KH) instabilities can affect this con- 
figuration. It has been shown that S PH techniques have 
difficulty handling KH instabilities (| Agertz et all I2007T ) 



where p c and p^ are the densities of the heavier and more 
dilute components, respectively, A is the spatial scale of 
the instability, and U is the relative velocity of the two 
phases. For a spherical cloud moving in a dilute medium 
the disruptio n time is 2 — 5 x tKH ( self-gravity can make it 
even longer) (jMurrav et al.l [19931 ). At the temperatures 
typical of penetrating filaments in the z ~ 0.5 MWH, 
T ~ lO 5 -^ which for approximate pressure equilibrium 
implies a density contrast of ~ 5 — 10. The overdense 
region in Figure [3] moves with velocities ~ 100 km/s 
through the hot medium. The largest available mode is 
responsible for the destruction of structures. The length 
of the structure is ~ 100 kpc and its typical width is ~ 20 
kpc. Assuming a 3 x tKH case for a wavelength compa- 
rable to the filament's width, it takes about 2 x 10 9 yrs 
for the instability to act. This is longer than the time 
over which the filament is destroyed by shocks and other 
instabilities and the time needed for gas to flow through 
the instability region visible in Figure [3J It is also pos- 
sible that the filament length is relevant which increases 
this timescale by factor of few. We therefore conclude 
that this structure can indeed penetrate into a hot halo 
to create such a configuration. Only much smaller modes 
can affect the infalling cold gas in a filament, perhaps en- 
hancing the surface shearing layer. 

Next, we estimat e the timescale for RT instabilities 
(jMurrav et al.lll993ft to act: 
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Fig. 3. — Example of the instability developing from an overdense filament in a hot dilute halo gas leading to cloud formation. Left panels 
show the projected density and the right panels show projected mass-weighted temperature in a box 200/i kpc on a side and lOO/i -1 kpc 
thick. 



where g is the gravitational acceleration. Taking the ac- 
celeration at 50 kpc in MWH and a density contrast ~ 5, 
we estimate the timescale for the instability at ~ 5 kpc 
scales, which corresponds to the example in Figure [3J 
t RT ~ 3.5 x 10 8 yr, comparable to t coo i ~ 3 — 4 x 10 8 
yr of the photo-ionized gas in a warm filament at T ~ 
1 - 1.2 x lO 5 ^ and n H = 2 x 10~ 4 cm~ 3 . Therefore, 
RT and cooling instabilities are acting together to form 
pressure supported clouds with T c ~ 1Q 4 K. Both of 
these timescales are shorter than the time over which 
the overdense structure is present in the inner halo. This 
process often simultaneously forms many clouds in the 
same region, yielding cloud complexes. At this stage, 
KH instabilities of similar wavelength could play a role. 
When clouds start to form, shearing velocities are slowed 
down by pressure gradients and shocks before they ac- 
celerate to terminal velocity so the role played by KH 
instabilities is unclear. 

4. DISCUSSION 

Halos formed in fully cosmological environments are 
complex and dynamic systems. In addition to the pro- 
cesses described above, infalling satellites can be stripped 
of gas as they pass through overdense filaments in the 
outskirts of halos, providing material for further insta- 
bilities and cloud formation. This effect can alleviate 
the need for a high density hot medium to produce large 
gaseous tails s uch as the one obser ved in the Magellanic 
Cloud system (jPutman et al.ll2003l ). An example of this 
is visible in the upper left corners of the bottom two 
panels in Figure [3J Infalling substructures create density 
waves and shocks in the halo, providing a source of fur- 
ther compression and acceleration which enables RT in- 
stabilities to develop even without the gravitational pull 
of the parent halo. Similarly, we expect density inver- 
sions and RT instabilities to be more common in sim- 



ulations with galactic outflows (both stellar and AGN), 
where winds can temporarily push, comp ress and acceler- 
ate b oth galactic and infalling material (|Di Matteo et al.l 
[2Q0l . 

The MWH has a growth history typical for its mass 
indicating that this cloud formation process is general. 
In future work we will discuss cloud formation in larger 
hal os where clouds fo rm at earlier times (e.g., see Figure 
6 in lKeres et al.|[2009l) . Halo clouds are likely part of halo 
absorption systems observed at intermediate redshift and 
some are likely analogs of the HVCs. However, from 
Figure Q] it is already clear that many of the cold clouds 
have velocities comparable to the galactic rot ation which 
will m ake searches for such objects difficult (jPeek et al.1 
2008). Cloud infall at large radii will also likely trigger 
star formation and provide kinematic enhancements in 
the outer disks. 

Cloud formation from cooling instabilities in an ide - 
alized setting was simulated in iKaufmann et all ((2009) , 
demonstrating that large entropy cores are needed for 
this mechani s m to operate. This is supported by 
iBinnev et al.l (|2009f) who examined relevant timescales 
for the onset of cooling instabilities that begin with small 
density enhancements and concluded that unless the en- 
tropy profile of the halo gas is shallow, buoyancy and 
conduction will prevent cloud formation. In our simu- 
lation clouds are seeded at intermediate radii starting 
from already moderate overdensities, for which this lin- 
ear analysis does not apply. However, the inner gas en- 
tropy and density profiles of MWH are quite shallow, so 
it is not clear if slower cloud formation is prevented. In- 
terestingly, clouds just below our resolution limit, < 32 
particles, do have a larger contribution from a hot gas 
component. More work is needed to determine if this is 
a physical effect or purely numerical if unresolved objects 
partially mix in hotter gas. In reality, even a moderate 
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metallicity of halo gas could increase cloud formation ef- 
ficiency. 

Both SPH and Eulerian codes have shortcomings in 
following the instabilities relevant when clouds move in 
a hot m edium. New tec hniques, such as the moving mesh 
code of lSpringell <|2009), may avoid these shortcomings, 
providing new insights into the formation and evolution 
of HVCs. 

4.1. Survivability of clouds 

It is interesting to determine if clouds formed in MWH 
can survive their trip through the inner halo and rain 
onto the galaxy, providing fresh fuel for star forma- 
tion. In the absence of magnetic fields, thermal con- 
duction will directly evaporate clouds smaller than about 
M c ~ 1O 5 M in a free fall time for the typical radii where 
clouds form in MWH (|Mo fc Miralda-Escudd[l99l . The 
evolution of more massive clouds is likely governed by 
a combination of RT and KH instabil ities (e.g. see 
iMurrav et ail 119931 : iMurrav fc Lm] 12004 . Our simula- 
tions do not have sufficient resolution to follow instabil- 
ities on such small scales and SPH itself might not be 
suitable for this task. 

We can, however, estimate relevant timescales. Ide- 
alized simulations show that clouds frag ment when they 
sweep up a mass comparable to their own (jMurrav fe Lin! 
12001 . Therefore, we integrate the descent of clouds on 
radial orbits, assuming they form at 70 kpc and need to 
reach 5 kpc to fall onto a galaxy. We use T/j = 10 6 K, 
T c = 1Q A K and a density profile p oc r^ 13 . We take 
(3 = 0.65 — 0.9 (bracketing the inner density profile of 
the MWH; for simplicity we assume Th = const.), and 
normalize the density at 70 kpc. 

We find that clouds above 2 - 3 x 10 5 M Q survive if 
clouds are compressed during their descent. This mass 
is ~ 5 times higher if the cloud radius remains fixed dur- 



ing infall. In both cases this is smaller than masses of 
the clouds formed in our simulation. We also directly 
integrate cloud orbits accounting for gravity and ram 
pressure, assuming that clouds have angular momentum 
similar to the disk gas. This only mildly increases the 
mass limit for cloud survival compared to radial infall. 
Therefore, massive clouds seeded in the inner halo are 
likely able to bring cold material to a central disk. In 
the non-wind simulations at z=0 this yields 0.6M Q /yr of 
fresh gas accretion, but if galactic outflows are included 
this increases to ~ IMq/yt (because more material is 
saved for later times). For the MW, whose mass is factor 
of ~ 2 larger than the simulated MWH, the rates are 
likely higher, providing enough fuel for continuous star 
formation over a several Gyr period. 

5. CONCLUSIONS 

Using a state-of-the-art numerical simulation of struc- 
ture formation in a ACDM universe we follow the for- 
mation of small ~ 10 6 Af Q cold, pressure confined clouds 
in a halo similar to the MW. Most clouds are not ran- 
domly formed from small fluctuations in a hot, T ~ 10 6 , 
medium as has often been supposed. Instead, they orig- 
inate from overdense, T ~ 10 5 , filamentary gas within a 
hot halo at low redshift. Such gas creates a density and 
entropy inversion in a gravitational field, driving simul- 
taneous cooling and RT instabilities. During their initial 
evolution, many clouds form from cometary structures 
while later evolution cannot be followed properly with 
our computational scheme. However, simple estimates 
indicate that a large fraction of such clouds could indeed 
survive their trip to the halo center and provide a supply 
of gas for late time star formation in galaxies. 

We are grateful to V. Springel for allowing us to use 
his code and to G. Besla, Y. Birnboim, A. Burkert and 
T. Cox for useful comments. 
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